____ _ _ _ _
| _ \ ___ | |_ (_) _ __ ___ __| | (_) __ _
| |_) | / _ \ | __| | | | '_ \ / _ \ / _| | | | / _ |
| _ < | __/ | |_ | | | |_) | | __/ | (_| | | | | (_| |
|_| \_\ \___| \__| |_| | .__/ \___| \__,_| |_| \__,_|
|_|
- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b- `b
Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―Β―
Backward-Algorithmus
ββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββ
top
Der Backward-Algorithmus (auch RΓΌckwΓ€rts-Algorithmus, RΓΌckwΓ€rts-Prozedur) berechnet mit Hilfe von Backward-Variablen die Wahrscheinlichkeit, in einem gegebenen Hidden-Markov-Modell (HMM) eine bestimmte Symbolsequenz zu beobachten. Der Algorithmus verwendet die Programmiermethode der dynamischen Programmierung.
Contents
β’ Markov-Modell
β’ Algorithmus
β’ KomplexitΓ€t
β’ Siehe auch
β’ Literatur
β’ Weblinks
ββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββ
Markov-Modell
Gegeben sei ein HMM Ξ» Ξ» = ( S ; V ; A ; B ; Ο Ο ) {\displaystyle \lambda =(S;V;A;B;\pi )} , wobei
β’ S {\displaystyle S} die Menge der verborgenen ZustΓ€nde,
β’ V {\displaystyle V} das Alphabet der beobachtbaren Symbole,
β’ A {\displaystyle A} die Γbergangsmatrix,
β’ B {\displaystyle B} die Matrix der Emissionswahrscheinlichkeiten,
β’ Ο Ο {\displaystyle \pi } die Anfangswahrscheinlichkeitsverteilung fΓΌr die mΓΆglichen AnfangszustΓ€nde,
bezeichnet.
Aufgabenstellung und Backward-Variablen
Gegeben sei ein Wort o = o 1 o 2 β¦ β¦ o T β β V β β {\displaystyle {\boldsymbol {o}}=o_{1}o_{2}\dots o_{T}\in V^{*}} . Der Backward-Algorithmus berechnet nun P ( o | Ξ» Ξ» ) {\displaystyle P({\boldsymbol {o}}|\lambda )} , also die Wahrscheinlichkeit, im vorhandenen Modell Ξ» Ξ» {\displaystyle \lambda } tatsΓ€chlich die Beobachtung o {\displaystyle {\boldsymbol {o}}} zu machen.
DafΓΌr werden die Backward-Variablen Ξ² Ξ² t ( i ) {\displaystyle \beta _{t}(i)} verwendet, sie bezeichnen die Wahrscheinlichkeit, das Suffix o t + 1 o t + 2 β¦ β¦ o T {\displaystyle o_{t+1}o_{t+2}\ldots o_{T}} zu beobachten, falls das HMM zum Zeitpunkt 1 β€ β€ t β€ β€ T {\displaystyle 1\leq t\leq T} im Zustand s i β β S {\displaystyle s_{i}\in S} gewesen ist:
Ξ² Ξ² t ( i ) = P ( o t + 1 o t + 2 β¦ β¦ o T | q t = s i ; Ξ» Ξ» ) {\displaystyle \beta _{t}(i)=P(o_{t+1}o_{t+2}\dotsc o_{T}|q_{t}=s_{i};\lambda )}
Algorithmus
Die Backward-Variablen werden rekursiv bestimmt:
Initialisierung
Ξ² Ξ² T ( i ) = 1 , 1 β€ β€ i β€ β€ | S | {\displaystyle \beta _{T}(i)=1,\qquad 1\leq i\leq \left|S\right|}
Rekursion
Ξ² Ξ² t ( i ) = β β j = 1 | S | b j ( o t + 1 ) β
β
a i j β
β
Ξ² Ξ² t + 1 ( j ) , 1 β€ β€ i β€ β€ | S | , 1 β€ β€ t < T {\displaystyle \beta _{t}(i)=\sum _{j=1}^{\left|S\right|}b_{j}(o_{t+1})\cdot a_{ij}\cdot \beta _{t+1}(j),\qquad 1\leq i\leq \left|S\right|,\ 1\leq t<T}
Termination
P ( o | Ξ» Ξ» ) = β β j = 1 | S | Ο Ο j β
β
b j ( o 1 ) β
β
Ξ² Ξ² 1 ( j ) {\displaystyle P({\boldsymbol {o}}|\lambda )=\sum _{j=1}^{\left|S\right|}\pi _{j}\cdot b_{j}(o_{1})\cdot \beta _{1}(j)}
KomplexitΓ€t
Die Matrix aller Backward-Variablen braucht O ( | S | β
β
T ) {\displaystyle O(|S|\cdot T)} Speicher, werden die Zwischenergebnisse im Anschluss nicht mehr verwendet, so reduziert sich der Platzbedarf auf O ( | S | ) {\displaystyle O(|S|)} , da nur mehr zwei Spalten der LΓ€nge | S | {\displaystyle |S|} benΓΆtigt werden, um die Werte von Ξ² Ξ² t + 1 ( i ) {\displaystyle \beta _{t+1}(i)} und Ξ² Ξ² t ( i ) {\displaystyle \beta _{t}(i)} in jedem Rekursionsschritt zu speichern.
FΓΌr jede einzelne Variable wird ΓΌber | S | {\displaystyle |S|} Zeilen summiert, also liegt die Laufzeit in O ( | S | 2 β
β
T ) {\displaystyle O(|S|^{2}\cdot T)} .
Weitere Anwendungen
Die Backward-Variablen Ξ² Ξ² t ( i ) {\displaystyle \beta _{t}(i)} werden zusammen mit den Forward-Variablen Ξ± Ξ± t ( i ) = P ( o 1 , o 2 , β¦ β¦ , o t , q t = s i | Ξ» Ξ» ) {\displaystyle \alpha _{t}(i)=P(o_{1},o_{2},\ldots ,o_{t},q_{t}=s_{i}|\lambda )} fΓΌr den Baum-Welch-Algorithmus zur LΓΆsung des mit Hidden-Markov-Modellen gegebenen Lernproblems benΓΆtigt.
AuΓerdem ermΓΆglicht deren Kenntnis die Bestimmung der Wahrscheinlichkeit bei der Beobachtung von o {\displaystyle {\boldsymbol {o}}} zu einem festen Zeitpunkt t {\displaystyle t} im Zustand s i {\displaystyle s_{i}} gewesen zu sein, denn nach dem Satz von Bayes gilt:
P ( q t = s i | o ; Ξ» Ξ» ) = Ξ± Ξ± t ( i ) β
β
Ξ² Ξ² t ( i ) P ( o | Ξ» Ξ» ) {\displaystyle P(q_{t}=s_{i}|{\boldsymbol {o}};\lambda )={\frac {\alpha _{t}(i)\cdot \beta _{t}(i)}{P({\boldsymbol {o}}|\lambda )}}}
Siehe auch
Literatur
β’ Richard Durbin, Sean R. Eddy, Anders Krogh, Graeme Mitchison: Biological sequence analysis. Probabilistic models of proteins and nucleic acids. 11th printing, corrected 10th reprinting. Cambridge University Press, Cambridge u. a. 2006, ISBN 0-521-62971-3, S. 59β60.
Weblinks
β’ Ernst G. Schukat-Talamazzini: Spezielle Musteranalysesysteme. (PDF, 1,3 MB). Vorlesung im Wintersemester 2012/2013 an der UniversitΓ€t Jena. Kapitel 5, Folie 34 ff.